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ABSTRACT 

We analyze an X-ray flux-limited, complete sample of 93 AGN at z<0.1, selected from the ROSAT 
Bright Survey. Two thirds of the sample are Seyfert 1 galaxies (Syl) and one third are Narrow-Line 
Seyfert 1 galaxies (NLSyl). We have obtained optical images of all objects. By modehng the host 
galaxy and the AGN central component we decompose the optical emission into nuclear, bulge and 
disk components, respectively. We find that the nuclear optical luminosity, thought to be associated 
with the accretion disk surrounding the active black hole, correlates with the X-ray luminosity, the 
radio luminosity and the black hole mass. 

Subject headings: galaxies: Seyfert — galaxies: nuclei — x-rays: galaxies 



1. INTRODUCTION 

The last few years have brought fundamental progress 
in understanding the suspected connection between ac- 
tive galactic nuclei (AGN) , characterized by non-thermal 
emission from a massive (10^ — 10^ M©) black hole (BH), 
and normal galaxies, whose extended emission is ther- 
mal. In particular, now 1) there is robust evidence that 
quasars reside in galaxies (Bahcall et al. 1997); 2) it is es- 
tablished that massive "quiescent" BHs are also present 
in the central regions of normal galaxies (Magorrian et 
al. 1998); 3) a statistical correlation between BH mass 
and mass or velocity dispersion of the host galaxy's bulge 
exists (Gebhardt et al. 2000; Ferrarese & Merritt 2000), 
which suggests a common evolution for galaxies and their 
central, massive BHs; 4) the previous correlations hold 
for normal as well as AGN-host galaxies (McLure & Dun- 
lop 2002). 

For AGN with typical masses of 10^ - 10^ Mq , the bulk 
of the accretion radiation is emitted in the soft X-ray and 
EUV/UV range. A fraction of this optically thick radi- 
ation is Compton-upscattered by a corona of electrons 
surrounding the BH. This process is thought to produce 
the typical power law X-ray spectra observed in all AGN. 
The primary soft X-ray/EUV radiation of the accretion 
disk also ionizes the surrounding gas, producing the well- 
known AGN emission lines. At radio frequencies, the 
core luminosity of the AGN is due to synchrotron emis- 
sion by relativistic electrons immersed in a seemingly 
global magnetic field. It is conceivable that these elec- 
trons are the same which Compton-upscatter the disk 
photons. It is therefore interesting to ask whether there 
is a direct correlation between radio, optical and X-ray 
emission in the central part of the AGN where the accre- 
tion disk is located. 

While information in the radio and X-ray bands is 
already available from all-sky surveys like the NRAO 
VLA Sky Survey (NVSS, Condon etal. 1998), FIRST 
(Faint Images of the Radio Sky at Twenty centimeters; 
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Becker et al. 1995) and the ROSAT Bright Survey (RBS, 
Schwope etal. 2000), optical data for the central kilo- 
parsec of AGN are more difficult to obtain. Indeed, in 
modeling AGN and their host galaxies, the information 
contained in the central pixels has been usually rejected 
(see e.g. Virani etal. 2000). With a new well-tested, 
well-performing and robust fitting algorithm (Kuhlbrodt 
etal. 2001), it has recently been possible to model both 
the host galaxy and the nucleus from optical images. 

Here we apply this decomposition algorithm to a new 
sample of objects, and use the results of the decomposi- 
tion to study the relations between the properties of the 
central point source, in particular the X-ray (0.5 — 2.0 
keV), optical (R band) and radio (1.49 GHz) emission. 
After presenting the sample selection in §1 and the cat- 
alog data in §2, the relations between the different pa- 
rameters are described in §3, and discussed in §4. Unless 
specified differently, a value of Ho=50 km Mpc^^s^"'^ has 
been used throughout this paper. 

2. CONSTRUCTION OF THE SAMPLE 

We used the ROSAT Bright Survey (RBS, Schwope 
etal. 2000), consisting of the ^ 2000 brightest, opti- 
cally identified sources detected during the ROSAT All 
Sky Survey (RASS, Voges etal. 1999). Selecting at X- 
rays has the advantage that the contribution of the emis- 
sion from the host galaxy is small with respect to the 
AGN. We selected the strongest (ROSAT PSPC count 
rate (cr) above 0.3 cts s~^ in the 0.5-2.0 keV range), local 
(0.009<z<0.1) AGN in order to achieve reasonable spa- 
tial resolution with ground-based optical observations. 
The sample consists of 93 AGN (black symbols in Fig. 

At the time when the work started, all of them were 
classified in NED^ as Syl galaxies, but in the meantime 
about 30% of them have been found to be NLSyl (Veron- 
Cetty & Veron 2001). The size of the sample thus allows 
us to work with 2 independent, complete samples, Syl 
(64) and NLSyl (29), with comparable X-ray luminosity. 
In both cases the observer has a direct view of the AGN, 
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i.e. unobstructed by the torus. 
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Fig. 1. — AGN in the RBS catalogue. Solid lines define the 
selection criteria in redshift (0.009<z<0.1; vertical lines) and in 
flux (cr > 0.3 cts/s; diagonal). Black symbols form the sample 
used in this work. 



3. DATA COLLECTION 

All optical data was acquired in dedicated observing 
runs, described in detail elsewhere (Salvato 2002; Salvato 
2003, in preparation). Here we give only some brief char- 
acteristics. In order to explore the relationship of nuclear 
luminosities in various wavelengths we collected data 
from public catalogs for the X-ray and radio wavelengths. 

Optical data. The entire sample has been observed 
in the R band with a typical exposure time of 10 
minutes at the 1.5m Danish telescope (Bessel R filter) 
at La Silla and of 20 minutes at the 1.2m telescope 
(Johnson R filter) at Calar Alto. We have applied the 
two-dimensional decomposition algorithm (Kuhlbrodt 
etal. 2001) to all images, which takes into account the 
point spread function (PSF) of the corresponding image. 
A three-component model has been applied, involving 
a point-like nuclear component, a bulge component 
with a R^/^ radial dependence, and an exponential disk 
component. Besides the spatial parameters we therefore 
obtain magnitudes of these three components. These 
derived magnitudes have been corrected for Galactic 
absorption, using the maps of Dickey & Lockman (1990) 
and converted to the visual extinction as derived by 
Predehl & Schmitt (1995). A comparison of our fit 
parameters (colors, disk and bulge parameters, etc.) 
with a sub-sample of galaxies studied already earlier 
(e.g. McLure & Dunlop 2001), reveals perfect agreement 
(details in Salvato 2002), and thus provides some 
confidence that also the nuclear light component is fitted 
correctly. In this optical nuclear luminosity is sup- 
posed to represent the accretion disk surrounding the BH. 

Black hole masses. BH masses are computed from 
the derived bulge luminosity, using a unique equation 
for Syl and NLSyl as in McLure & Dunlop (2002): 
;o5(Mbh/Mq)= -0.50(±0.02)Mr - 2.96(±0.48) which 



is valid under the assumption that differences between 
Syl and NLSyl are due to geometrical rather than 
physical effects. Our masses (Tab. QJ, after correcting 
for the different cosmology, agree well with those derived 
by reverberation mapping (Kaspi et al. 2000) for the 
common sources. 

X-ray luminosity. Typically, non-active galaxies have 
very low X-ray luminosities; thus the X-ray emission 
at a ROSAT/PSPC count rate above 0.3 cts/s is likely 
produced only in the nuclei of AGN. However, the X-ray 
emission is not directly related to the accretion disk, 
because due to the canonical temperatures of the disk 
emission peaks in the UV. The X-ray luminosity has 
been computed using the X-ray flux in the 0.5-2.0 keV 
(Tab. 2 in Schwope etal. 2000), assuming a power law 
spectrum with photon index of -2 and applying the 
corresponding Galactic absorption and k correction. 
We caution that this approach yields unreliable and 
inconsistent results for few NLSyl in our sample, having 
very steep, soft X-ray spectra. 

Radio data. We cross-correlated our sample with the 
NVSS and FIRST radio catalogues, and, owing to their 
limited sky coverage, 53 matches were found. Both 
surveys are centered at 1.49 GHz (A = 20 cm), where 
the synchrotron emission of our targets may be as- 
sociated with the nuclear activity but also with the 
spatially-extended star formation activity, if any. NVSS 
and FIRST map the sky at different spatial resolution 
(45"and 5", respectively), and thus allow to disentan- 
gle the two potential contributions. In fact, if the unre- 
solved AGN produces all the radio emission, the fluxes 
listed in the two catalogues must be consistent; otherwise 
the NVSS flux must be larger. A flux-flux plot shows 
that galaxies with and without the extra contribution to 
the radio emission from star formation activity can be 
separated by a line of equation log(Li. 49GHz (NVSS) = 
log(Li.49GHz (FIRST)) -I- 0.3. This selection criterion is 
robust against a maximal 20% variability of the AGN 
radio flux (Falcke etal. 2001). Thus the radio luminosi- 
ties of 21 out of 53 galaxies are attributed only to their 
nuclear activity. 

Table 1 Comparison of BH masses obtained 

FOLLOWING McLuRE&DuNLOP (2002) AND VIA 

reverberation mapping (Kaspi etal. 2000). All 
values are computed using ho = 75 km mpc"^s"^. 



Source Name 


Kaspi et al. 


This work 




log{Mth/MQ) 


log(Mi,;,/M0) 


nOll (Fairall 9) 


7.59^ 8.01 


7.89^8.83 


n054 (PG 1211-^14) 


7.45 ^ 7.7 


8.44^ 9.34 


n065 (IC 4329A) 


7.25 i 


6.76 ^7.78 


n069 (NGC 5548) 


8.02 ^8.16 


7.14^ 8.12 



4. CORRELATION ANALYSIS 

Optical vs. X-ray emission. A strong linear relation is 
found between the optical accretion disk and the X-ray 
luminosity (Fig. 13), of the form 

\og{Lx)^a-\og{LR{Nucj)+b (1) 
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Table 2 Coefficients for the equations 1 and 2. The 
Spearman Rank correlation coefficients p , are 

GIVEN IN COL. 3 and 6. 
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Equation 1 




Equation 2 




a b 


P 


c d 


P 


All 0.58±0.06 18.44±2.6 


0.74 


0.96±0.24 -4.18±10.7 


0.68 


Syl 0.56±0.07 19.48±3.0 


0.72 


0.74±0.22 4.80±9.6 


0.60 


NLSyl 0.61±0.11 17.09±4.9 


0.82 


1.18±0.52 -13.95±22.9 


0.92 



where the coefficients "a" and "b" have been deter- 
mined separately for the whole sample (All) and for 
the Syl and NLSyl galaxies alone, respectively (Tab. 
121). The consistency of the relations is confirmed by 
the Spearman Rank correlation coefficient, p, which 
is higher than 0.7. From the plot and the table it is 
obvious that Syl and NLSyl behave in the same way. 
The scatter is about ±0.5 dex. There are (at least) 
two obvious reasons for this scatter: first, optical and 
X-ray fluxes have not been measured simultaneously, 
and due to the known variability of Seyferts a factor of 
a few can easily be accounted for. Second, the X-ray 
fluxes have been computed on the basis of the total 
measured counts and an assumed power law slope. It 
would have been more correct (though not possible 
in most cases due to photon statistics) to fit also the 
power law slope, and thus obtain a more accurate 
X-ray flux. One may think that the Ha line may 
produce a non-negligible contribution to the continuum 
emission in the R-band. For testing this hypothesis, the 
template spectrum of a Syl has been convolved with the 
transmission curve of the R-band filter. We computed 
the dilution factor from Ha, for different [NII]/Ha ratios 
and for different EW(Hq,), also taking into account that 
Ha may be enhanced by galaxy interaction. The dilu- 
tion factor is not more than 1% and thus being negligible. 

Radio vs. X-ray emission. We also find that the X- 
ray emission correlates well with the radio continuum 
emission. As in the previous case, the relation is linear, 
and has been fit in the form 



fog(i 



FIRST 



(Nuc)) 



\og{Lx) + d 



(2) 



of Tab. 121 Here we used only those 21 galaxies for which 
we have information about nuclear radio emission. As 
previously, the linear regression has been computed for 
all objects, and for the Syl (14 galaxies) and NLSyl (7 
galaxies) independently. Due to the small size of the two 
subsamples the errors are quite high and the Spearman 
Rank correlation coefficient in not really reliable. Never- 
theless, it appears that NLSyl may have a steeper slope 
than the Syl galaxies. 

Optical vs BH masses. The optical luminosity of the 
accretion disk is also correlating with the BH mass com- 
puted according to §3.2 (Fig. The solid line indicates 
the predicted relationship for an optically thick accre- 
tion disk around a non-rotating BH (Shakura & Sunyaev 
1973), where for a given BH mass the R band flux de- 
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Fig. 2. — X-ray versus optical (R-band) luminosity for the sam- 
ple. Asterisks (open triangles) represent Syl (NLSyl) galaxies. For 
30 out of the 93 objects imaging was obtained under bad seeing 
condition (>2"). These are marked by open circles around their 
respective symbols and were ignored for determining the linear cor- 
relation between the X-ray and optical luminosities. Larger seeing 
implies more flux in the nuclear component relative to the bulge 
component, and therefore those encircled data points have to be 
treated as upper limits on the nuclear optical emission. At bottom 
left the expected error for the nuclear, optical luminosity is given. 
Syl and NLSyl behave in the same way. 
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Fig. 3. — Radio versus X-ray luminosity for the sample. Symbols 
are the same as in Fig. |^ AGN with lower count rate than those 
we considered (grey asterisks) are also plotted, with the intention 
to show that our relation is not driven by the Malmquist bias. 



pends on the temperature (and thus accretion rate M) of 
the object {M ~ O.lMEdd chosen here). The curvature of 
this line is basically a measure of the differing bolometric 
correction for varying BH mass (and thus temperature) . 

5. CONCLUSIONS 

We find that the optical luminosity of the nuclei of 
Seyfert galaxies correlates with the power-law X-ray 
emission, the radio emission from the core as well as the 
Black Hole mass. 

These relations are not driven by the Malmquist bias, 
which is introduced whenever a sample is selected via a 
fiux- limit. In fact, the additional 9 Syl galaxies with cr 
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log(MBH) [M,,„] 

Fig. 4. — R-band luminosity versus BH mass for the sample. 
Symbols are the same as in Fig. |21 At bottom left typical errors are 
shown. The plotted line is not a fit, but the predicted relationship 
for an optically thick accretion disk with M = 0.1 MEdd- 



< 0.3 cts/s and NVSS/FIRST counterpart confirm the 
general trend (see Fig. O grey asterisks.) 

The correlation of the optical nuclear luminosity with 
the X-ray luminosity could potentially be understood if 
the canonical blackbody-like continuum emission from 
the accretion disk is enhanced due to irradiation by soft 
X-rays, so that an increase in the X-ray flux would create 
a larger optical flux. This is a well-known effect in X-ray 
binaries, though it has not been possible so far to quan- 
tify this relation, as it depends on too many parameters 
(e.g. X-ray spectrum of the irradiating source, temper- 
ature and ionization state of the disk, geometry of the 
disk). 

Kriss etal. (1982) and Koratkar etal. (1995) found 
a strong correlation between the soft X-ray luminosity 
and the broad Hp and Ha line luminosity, respectively, 
in quasars and Syl galaxies. It has been argued that 
an increased accretion disc flux would not only cause a 
larger photoionization (and thus Hq flux), but also pro- 
vide more seed photons which would be up-scattered to 
produce the power law X-ray spectrum, and thus the X- 
ray flux. This effect is unlikely to explain our correlation, 
since simulations show that the contribution of the Hq 
emission line to the R band flux is less than 1%. 

Including the correlation with the radio luminosity 
leads to the consideration of a less direct, yet at least as 



plausible, coupling between the luminosities in the differ- 
ent wavelengths, namely via the electron corona which 
is supposed to surround the accretion disk. The elec- 
tron temperature distribution is expected to depend on 
the accretion disk properties, primarily the accretion rate 
and disk temperature. On the other hand, the tempera- 
ture and density of the corona determine the efficiency of 
the Compton scattering, and thus directly determine the 
slope and power of the X-ray spectrum. Thus, a change 
in e.g. the accretion rate would not only lead to a vary- 
ing optical emission, but via the corona also to a varying 
X-ray and radio emission. The physical relation between 
the processes which produce the emission in the X-ray 
and radio remains at the level of speculation. In mi- 
croquasars, the stellar-binary analogs of jet emitting ob- 
jects in our Galaxy, Rau & Greiner (2002) have recently 
found a similar correlation between the strength of the 
radio emission and the X-ray power law slope. Owing 
to the smaller timescales involved in microquasars (due 
to their smaller black hole mass), a whole continuum of 
jet strengths has been observed from bright, spatially re- 
solved ejections over "baby-jets" down to quasi-steady 
radio (synchrotron) emission. Adopting this similarity, 
one could argue that the core radio emission in Seyfert 
galaxies witnesses matter outflow from the central black 
hole. Objects with higher activity at X-ray wavelengths, 
i.e. higher accretion rate, may thus also show higher mass 
outflows from their central engine, thus implying a larger 
radio flux. Indeed, this is what has been found when 
comparing the X-ray and radio fluxes of galactic black 
hole binaries in the low/hard state: the radio flux scales 
as X-ray flux to the power of 0.7 (Gallo etal. 2003). 
Surprisingly, this is (within the errors) the same corre- 
lation as we find (coefficient c in Tab. |2l. Certainly, 
further investigations are necessary to understand the 
processes involved and their interrelation, but it is en- 
couraging that seemingly the same processes are at work 
in AGN and galactic black hole binaries. 
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